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ABSTRACT 

A  method  for  reconstructing  the  3D  shape  of  a  single  grain  in  a  polycrystal  from  far-held 
diffraction  data  is  presented.  The  reconstruction  is  performed  using  an  iterative  algorithm, 
algebraic  reconstruction  technique  (ART),  to  solve  a  linear  system  of  equations  obtained 
from  data  from  the  Boxscan  technique.  The  Boxscan  technique  is  a  3D  X-ray  diffraction- 
type  scanning  technique,  which  -  in  geometric  tomography  terms  -  provides  2-dimensional 
X-rays  of  the  object  as  data.  In  this  paper,  we  report  on  the  first  experimental  3D  grain 
shape  reconstruction  based  on  a  beta-titanium  sample.  By  comparing  our  results  to  high- 
resolution  phase  contrast  tomography,  we  find  that  the  average  error  on  the  position  of  the 
grain  boundary  in  this  sample  is  2.7  micrometers  for  a  grain  radius  of  17  micrometers.  This 
result  is  similar  to  other  current  methods. 


1.  INTRODUCTION 

The  last  decades  development  of  nondestructive  3D  characterization  techniques  for  poly¬ 
crystals  has  been  important  in  the  held  of  materials  science.  The  outcome  of  such  a  char¬ 
acterization  could  be  a  map  of  the  positions  of  the  grain  centers  in  3D  or  even  better:  a  full 
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3D  description  of  the  topology  in  the  volume  (which  is  usually  called  a  3D  grain  map). 
Within  the  general  approach  of  3-Dimensional  X-Ray  Diffraction  (3DXRD)  several  grain 
mapping  techniques  have  been  explored.  In  the  first  one  developed  (Lauridsen,  Schmidt, 
Suter  and  Poulsen  2001;  Poulsen  2004)  the  sample  is  illuminated  layer-by-layer  with  a  line- 
shaped  high-energy  X-ray  beam  and  the  diffracted  signal  is  recorded  on  both  a  near-field 
and  a  far-held  detector  as  the  sample  is  rotated.  Where  the  far-held  signal  gives  informa¬ 
tion  for  resolving  the  grain  orientation,  the  near-held  detector  shows  a  2D  projection  of  the 
diffracting  grains  in  the  illuminated  sample  layer.  As  shown  in  Poulsen  and  Fu  (2003)  and 
Fu,  Knudsen,  Poulsen,  Herman,  Carvalho  and  Liao  (2006)  applying  algebraic  reconstruction 
techniques  (ART)  to  the  near- field  signal  enables  the  possibility  to  produce  spacefilling  grain 
maps  of  the  illuminated  sample  layer. 

Ideas  from  3DXRD  and  absorption  tomography  imaging  triggered  the  development  of  Diffrac¬ 
tion  Contrast  Tomography  (DCT)  (Ludwig,  Schmidt,  Lauridsen  and  Poulsen  2008;  John¬ 
son,  King,  Honnicke,  Marrow  and  Ludwig  2008;  Ludwig,  Reischig,  King,  Herbig,  Lauridsen, 
Johnson,  Marrow  and  Buffiere  2009).  While  the  sample  is  illuminated  layer-by-layer  in  the 
methods  mentioned  above,  DCT  illuminates  the  entire  volume  of  interest  and  records  the 
absorption  and  the  diffraction  signal  on  the  same  detector  positioned  in  the  near-field  as  the 
sample  is  rotated.  Again,  the  detector  images  show  2D  projections  of  the  diffracting  grains 
for  different  rotation  angles.  By  applying  a  set  of  criteria  based  on  the  geometry  and  the 
diffraction  angles  it  is  possible  to  determine  the  2D  projections  associated  with  each  grain. 
The  3D  shape  of  the  grain  is  computed  from  the  2D  projections  using  ART.  The  use  of  data 
back  projection  is  common  to  the  3DXRD  and  DCT  making  both  methods  sensitive  to  spot 
overlap  on  the  detector. 

High  Energy  X-ray  Diffraction  Microscopy  (HEDM)  (Suter,  Hennessy,  Xiao  and  Lienert 
2006),  and  the  methods  presented  in  Alpers,  Poulsen,  Knudsen  and  Herman  (2006)  are  vari¬ 
ants  of  the  3DXRD  but  the  grain  maps  are  reconstructed  by  forward  simulating  candidate- 
solutions  in  a  Monte  Carlo  optimization  scheme.  The  problems  with  spot  overlap  are  over¬ 
come,  but  powerful  computer  clusters  might  be  necessary  for  the  methods  to  be  useful  in 
practice. 

The  spatial  resolutions  of  the  methods  discussed  above  are  dependent  on  the  pixel  size  of 
the  detector  and  the  methods  require  a  detector  placed  close  (1-10  mm)  to  the  sample.  Such 
detector  positioning  is  in  some  cases  prohibited  due  to  the  presence  of  spacious  experimental 
equipment,  e.g.  a  furnace  or  a  stress  rig.  In  this  case  only  a  far-field  detector  can  be  used. 
With  only  a  far-field  detector  in  the  set-up,  with  3DXRD  one  could  obtain  information 
on  the  crystallography,  stress-state,  volume  and  center  of  mass  of  more  than  1000  grains 
(Oddershede,  Schmidt,  Poulsen,  Sorensen,  Wright  and  Reimers  2010).  However,  the  spatial 
resolution  is  reduced  to  approximately  10  micrometers  and  no  shape  information  is  available. 
Yet  the  combination  of  position  and  volume  makes  it  possible  to  generate  approximate  grain 
maps  based  on  a  Laguerre  tesselation  (Lyckegaard,  Lauridsen,  Ludwig,  Fonda  and  Poulsen 
2010).  Laguerre  tesselations  are  extensions  of  the  Voronoi  tesselation  and  are  found  to  pro¬ 
vide  a  better  approximation  to  the  underlying  grain  structure. 

In  the  techniques  discussed  above,  the  limiting  factor  has  been  the  detector  resolution. 
Where  the  detectors  have  resolutions  in  the  range  of  micrometers,  the  scanning  mechanics 
has  sub-micrometer  resolution,  some  even  down  to  10s  of  nanometers.  Thus,  using  the  me¬ 
chanics  to  probe  the  sample  by  translating  vertically  and  horizontally  through  an  extended 
line  beam  while  recording  the  diffracted  X-rays  on  a  far-field  detector  could  overcome  the 
problem  of  limited  detector  resolution.  A  technique  based  on  this  principle  has  been  de¬ 
veloped  through  the  Boxscan  technique  (Lyckegaard,  Poulsen,  Ludwig,  Fonda,  Margulies, 
G5tz,  Sorensen,  Dey  and  Lauridsen  2010).  The  Boxscan  technique  is  a  variant  of  3DXRD, 
but  uses  only  a  far-field  detector.  It  is  a  scanning  technique  where  two  datasets  are  ac¬ 
quired:  one  where  the  sample  is  scanned  vertically  and  one  where  it  is  scanned  horizontally. 
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By  combining  the  two  datasets  and  using  the  knowledge  on  the  diffracting  angles  and  the 
geometry,  the  size  and  the  positions  of  the  grain  centers  in  the  sample  can  be  determined 
with  high  precision.  In  addition  to  this,  several  projections  (2D  X-rays  in  terms  of  geomet¬ 
ric  tomography)  of  the  individual  grains  are  obtained.  As  before,  the  position  of  the  grain 
center  and  the  size  enables  the  construction  of  a  grain  map  using  a  Laguerre  tesselation. 
In  this  paper  we  point  out  that  boxscan  data  comprise  additional  information,  namely  ID 
profiles  of  the  grain  shapes.  We  discuss  how  to  exploit  this  information  and  report  on  the 
first  successful  3D  grain  shape  reconstruction  directly  from  Boxscan  data  (without  the  use 
of  a  tesselation). 


2.  BOXSCAN  TECHNIQUE 

2.1  Experimental  setup.  The  experimental  setup  for  the  Boxscan  technique  is  sketched  in 
figure  1.  The  sample  is  mounted  on  top  of  a  rotation  stage  u  and  a  set  of  stages  for  trans¬ 
lations  in  x,y,z  (laboratory  coordinatesytem) .  u  rotates  around  the  2- axis.  The  incoming 
focused  X-ray  beam  is  dimensioned  by  the  motorized  slits  located  in  front  of  the  sample. 
The  diffracted  X-rays  are  recorded  by  a  2D  detector  located  in  the  far-field,  i.e.  no  spatial 
information  of  the  grain  shape  is  recorded. 

It  is  convenient  to  introduce  a  sample  coordinatesytem  ( xs,ys,zs )  which  is  fixed  to  the 
sample,  thus  its  rotation  with  respect  to  the  laboratory  coordinatesytem  ( x,y,z )  is  u>.  For 
u  =  0,  (xs.  ys ,  2.5 )  and  (x,  y,  z)  coincide.  Eventually,  the  grain  positions  will  be  expressed  in 
the  sample  coordinate  system. 


Fig.  1:  Experimental  setup  of  the  Boxscan  technique  with  definition  of  angles  20,  lu 
and  q. 


2.2  Data  acquisition.  The  data  acquisition  consists  of  two  parts,  one  where  the  sample  is 
translated  vertically  in  2  and  one  where  it  is  translated  horizontally  in  y.  In  the  first  case, 
the  beam  is  slit  down  vertically  to  a  horizontal  line  beam.  For  each  vertical  translation 
2  =  2i...2c,  a  set  of  detector  images  are  recorded  for  a  series  of  sample  rotations  u  =  u>i...uja. 
In  the  second  case,  the  beam  is  slit  down  horizontally  to  a  vertical  line  beam.  For  each 
horizontal  translation  y  =  y\...yb,  a  set  of  detector  images  are  recorded  for  the  same  series 
of  sample  rotations  c 0  =  uq.-.u oa. 

2.3  Data  analysis.  The  diffraction  spots  are  extracted  from  the  images  and  the  values  (u,  20,  //, 
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z)  or  (u,  26,  rj,  y )  of  the  diffraction  angle  and  the  vertical  or  horizontal  translations  are  stored 
in  a  database.  The  subsets  of  diffraction  spots  with  similar  values  (lo,  26,  rj)  and  neighbouring 
values  in  y  or  z  are  denoted  a  profile,  see  figure  2.  Each  profile  corresponds  to  the  projection 
of  the  diffracting  grain’s  shape  onto  the  axis  of  translation  in  the  laboratory  coordinate  sys¬ 
tem.  Likewise,  the  intensity  weighted  average  of  the  profile  is  the  grains  centre  of  mass  on 
the  axis  of  translation  in  the  laboratory  coordinate  system,  ycms(co)  or  zcms(uj),  for  the  given 
upvalue.  It  should  be  stressed  here  that  the  grain  positions  are  fixed  in  the  sample  coordi¬ 
nate  system  but  moves  as  a  function  of  u  in  the  laboratory  coordinate  system.  Each  profile 
and  the  value  (u,  26,  rj,  ycms( &))  or  (a,  26,  rj,  zcms(uo ))  is  stored  in  the  database.  Combining  a 
horizontal  with  a  vertical  profile,  (u>,26,rj,ycms(u))  and  (u,  26,  rj,  zcms(uo)),  where  (u,29,rj) 
are  the  same,  i.e.  from  the  same  grain,  gives  a  set  of  matches  (u>,26,rj,ycrns(uj),  zcms(u>)) 
which  can  be  used  for  indexing  the  grains  in  the  dataset. 


Fig.  2:  Generation  of  vertical  profile  for  one  diffraction  spot.  As  the  sample  is  trans¬ 
lated  with  respect  to  the  beam  (right)  the  diffraction  spot  changes  appearance  (mid¬ 
dle).  The  profile  is  defined  as  the  integrated  intensity  in  the  spot  as  function  of  z 
(left). 


2.4  Indexing.  For  each  match,  the  value  (u>,  ycms(uj),  zcms(u))  maps  to  a  line  parallel  to  the 
(xs,  ys)-plane  in  the  3D  sample  coordinate  system  and  (cu,  26,  rj)  maps  to  a  line  in  3D  orien¬ 
tation  space,  represented  by  Rodrigues  vectors  (for  a  definition  of  Rodrigues  space  see  e.g. 
Morawiec  and  Field  (1996)).  A  subset  of  the  matches  which  fulfils  zcms(co )  =  z0  ±  A z  are 
extracted  from  the  database  and  the  intersection  points  (xs,  ys,  zs)  in  the  sample  coordinate 
system  and  (ri,r2,r3)  in  Rodrigues  space  for  all  possible  combination  of  two  matches  in 
the  subset  are  computed.  If  the  intersection  point  ( xs,ys,zs )  is  inside  the  scanned  volume 
and  both  exist  in  the  Rodrigues  space  and  is  inside  the  fundamental  zone,  it  is  accepted 
and  added  to  a  list  of  6D  points  (xs,  ys ,  zs  ,r\,  r2 ,  r3)  for  the  subset.  Now,  by  definition  each 
grains  centre  of  mass  and  crystallographic  orientation  is  described  by  a  single  point  in  the 
6D  space  {xs,  ys,  zs,  rq,  r2,  r3).  Two  matches  from  the  same  grain  will  intersect  in  exactly 
this  point,  thus  the  set  of  6D  points  for  the  subset  will  have  a  high  density  of  points  where 
a  grain  possibly  resides  in  6D.  To  find  grain  candidates  in  the  set  of  6D  points,  a  k- nearest 
neighbours  search  is  performed.  Points  in  6D  with  high  density  are  simulated  and  compared 
with  the  measured  data.  If  there  is  agreement  between  the  simulated  grain  candidates  and 
measured  data,  the  grain  candidate  is  accepted  as  a  grain. 

2.5  Verification.  The  Boxscan  technique  was  applied  to  a  cylindrical  beta  titanium  sample 
(300  micrometers  in  diameter)  with  an  average  grain  size  of  25  micrometers.  The  sample 
was  scanned  in  a  vertical  range  of  150  micrometers  with  a  step  size  of  5Z  =  4  micrometers 


3D  grain  reconstruction  from  boxscan  data 


and  in  a  horizontal  range  of  330  micrometers  with  a  step  size  of  Sy  =  5  micrometers.  The 
beam  had  a  vertical  size  of  20  micrometers  for  translations  in  z  and  a  horizontal  size  of 
30  micrometers  for  translations  in  y.  The  sample  was  probed  in  u  in  two  perpendicular  30 
degrees  intervals.  The  experimental  setup  used  prohibited  scanning  the  full  horizontal  width 
of  the  sample,  allowing  only  a  diamond  shaped  region  inside  the  sample  to  be  analysed. 

In  order  to  validate  the  grain  centers  determined  by  the  Boxscan  technique,  they  were 
compared  to  grain  centers  derived  from  a  phase  contrast  tomography  scan  (PCT)  of  the 
same  sample  (Cloetens,  Ludwig,  Baruchel,  Guigay,  Pernot-Rejmankova,  Salome-Pateyron, 
Schlenker,  Buffiere,  Maire  and  Peix  1999).  The  grains  in  the  imaged  volume  were  segmented 
by  a  series  of  image  processing  steps.  First  the  grain  boundaries  were  enhanced  and  noise 
suppressed  by  correlating  with  a  model  of  a  grain  boundary.  After  this  followed  a  bina- 
rization,  Euclidean  distance  transform  and  a  water  shedding.  After  the  binarization  a  few 
manual  corrections  were  needed  to  segment  the  grains  correctly. 

The  grain  centers  were  computed  for  the  grains  found  from  PCT  and  compared  to  the  ones 
from  the  Boxscan  technique.  The  Boxscan  technique  determined  133  of  159  verification 
grains  correctly.  The  grain  that  were  not  found  were  either  on  the  limit  of  the  scanned 
region  or  very  small  grains.  The  error  on  the  grains  centers  was  2.6  micrometers,  i.e.  the 
precision  of  the  grain  centers  were  better  than  the  step  sizes  Sz  =  4  micrometers  and  Sy  —  5 
micrometers. 


3.  APPLICATION  OF  ALGEBRAIC  RECONSTRUCTION  TECHNIQUE  TO 

BOXSCAN  DATA 

The  algebraic  reconstruction  technique  (ART)  as  described  by  Kak  and  Slaney  (1988)  is 
a  framework  for  solving  m  linear  equations  in  n  unknowns  Ax  =  b,  with  dimensions 
Amxn,  xnxi,  bmxi.  In  the  random  version  of  ART,  we  iteratively  produce  solutions  x*,+i 
by  projecting  the  current  solution  X/.  onto  a  random  row  of  the  A-matrix: 

b?  —  <pxt. 

Xfc+1  =  Xfc  +  A— - — rr^—SLi 

INI 

where  a,  correspond  to  the  ith  row  of  the  A-matrix  and  A  is  a  relaxation  parameter. 

The  Boxscan  technique  gives  shape  information  of  the  individual  grains  through  the  profiles. 
For  a  point  (t,  I )  in  a  profile,  where  I  is  the  measured  intensity  and  t  is  a  position  in  y  or  z,  I 
scales  linearly  with  the  integral  of  the  grain  volume  which  the  beam  illuminated  at  position 
t.  The  beam  positions  y  and  z  and  the  beam  widths  are  known,  hence  it  is  possible  to  set 
up  a  system  of  linear  equations  by  subdividing  the  sample  space  (x,  y,  z )  into  n  voxels.  The 
voxels  in  the  subdivided  sample  space  are  given  a  weight  between  0  and  1  according  to  the 
fraction  of  the  voxel  that  is  illuminated  by  the  beam  and  placed  into  the  A-matrix  as  a  row. 
Similarly,  /  is  stored  in  the  same  row  index  of  the  b- vector.  This  is  repeated  for  all  beam 
positions  in  all  profiles  belonging  to  one  grain.  The  unknown  x-vector  holds  the  intensity 
of  the  voxels  and  will  for  a  correct  solution  give  a  3D  description  of  the  grain. 

Assuming  that  the  resolution  is  such  that  we  can  assume  that  the  voxels  are  either  completely 
void  or  filled  by  a  grain,  it  is  natural  to  consider  x-vectors  that  are  either  0-  or  1-valued. 
ART,  however,  typically  returns  solutions  x  that  are  real  valued,  and  not  necessarily  binary. 
A  common  way  to  overcome  this  problem  is  to  introduce  a  threshold.  In  the  following,  we 
find  a  suitable  threshold  value  by  using  the  knowledge  obtainable  from  the  profiles.  Since 
each  profile  is  a  projection  of  the  grain  volume,  it  is  possible  to  estimate  the  grain  radius 
from  the  profiles  belonging  to  a  grain.  With  a  sphere  as  grain  model,  a  projection  of  the 
grain  convoluted  with  the  beam  profile  can  be  fitted  to  all  profiles  in  a  grain.  The  average 
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of  all  the  fitted  sphere  radii  is  used  for  calculating  the  grain  volume,  which  is  used  to  set 
the  threshold  on  the  ART  solution  such  that  the  volume  of  the  binarized  ART  solution  is 
equal  to  the  estimated  grain  volume. 


4.  RESULTS 

For  the  chosen  grain  in  the  Boxscan  dataset,  the  sample  space  was  voxelized  with  a  voxel  size 
of  2  micrometers.  This  in  conjuction  with  voxel  weights  and  profiles  was  used  to  generate 
the  linear  system  Ax  =  b.  The  system  had  15625  unknowns  and  931  equations  making  it  an 
under-determined  system.  The  ART  algorithm  was  run  with  A  =  10-5  and  terminated  when 
the  figure  of  merit  ||Ax  —  b||/||b||  had  converged.  The  average  grain  radius  was  estimated 
to  17  micrometers  and  the  corresponding  threshold  was  set  as  described  above.  In  figure  3 
the  reconstructed  binarized  grain  is  compared  to  the  same  grain  segmented  from  the  PCT 
verification  data  resampled  to  a  2  micrometers  resolution. 


Fig.  3:  Left:  3D  grain  shape  found  with  ART.  Right:  3D  grain  shape  from  verification 
data.  The  units  on  the  axes  are  micrometers. 


The  grain  reconstructed  with  ART  is  not  as  detailed  as  the  verification  data  but  the  overall 
shape  is  caught.  The  ART  solution  labels  87%  of  the  voxels  in  the  verification  data  correctly. 
For  a  better  comparison  of  the  ART  solution  to  the  verification  data,  the  distance  between 
the  grain  boundaries  of  the  two  grains  are  mapped  onto  the  surface  of  the  grain  from  the 
verification  data,  see  figure  4,  left.  As  the  figure  shows,  the  smooth  ART  solution  is  not 
capable  of  reconstructing  areas  around  corners.  The  average  distance  between  the  two  grain 
boundaries  is  2.7  micrometers.  Figure  4,  right,  shows  a  single  layer  of  the  two  grains  where 
red  is  verification  data,  green  Boxscan  data  and  yellow  is  overlap.  Again  it  is  possible  to  see 
that  the  ART  solution  has  a  smoother  and  more  rounded  shape  than  the  verification  data. 


5.  DISCUSSION 

As  laid  out  in  the  introduction,  the  question  to  be  answered  was  if  additional  projection 
information  could  improve  the  3D  grain  shape  reconstructions.  Based  on  the  Laguerre 
tessellation  approach,  in  Lyckegaard  et  al.  (2010a)  the  number  of  correctly  assigned  voxels 
was  measured  to  be  86%  for  a  sample  similar  to  the  one  used  here.  The  86%  was  measured 
as  an  average  over  all  grains  where  the  grain  neighbour  relationship  was  known.  In  our  case 
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Fig.  4:  Left:  Distance  between  the  grain  boundaries  of  the  ART  reconstruction  and 
the  verification  data  for  the  morphology  of  the  same  grain.  The  colour  scale  is  in  units 
of  micrometers.  Right:  A  superposition  of  the  same  layer  in  the  two  reconstructions 
of  the  grain.  Here  red  is  verification  data,  green  Boxscan  data  and  yellow  is  overlap. 
The  pixel  size  is  2  micrometers. 


where  the  ART  was  used  for  reconstructing  the  3D  grain  shape,  87%  of  the  grains  voxels 
were  assigned  correctly  compared  to  the  verification  data.  This  indicates  that  the  use  of 
the  additional  projection  information  gives  similar  results  even  without  information  about 
neighbouring  grains. 

Since  grains  are  positioned  in  a  spacefilling  network  and  thus  are  further  constrained,  it  may 
be  possible  to  improve  the  grain  shape  reconstruction  even  further  by  using  this  knowledge. 
In  other  words,  the  current  result  from  ART  reconstruction  may  be  improved  by  extending 
the  ART  implementation  to  simultaneously  solve  the  3D  grain  shape  of  a  number  of  grains, 
or  by  combining  the  output  of  the  ART  reconstruction  with  a  Monte  Carlo  scheme  where  the 
shape  of  a  number  of  grains  are  refined  simultaneously  by  a  forward  simulation.  Another 
and  maybe  more  obvious  way  to  improve  the  result  would  be  to  vary  the  experimental 
settings.  Extending  the  cj-range  such  that  more  projections  are  measured  or  decreasing  the 
step  sizes  Sy  and  6Z  could  also  improve  the  result. 


6.  CONCLUSIONS 


A  method  for  reconstructing  the  3D  shape  of  single  grains  from  Boxscan  data  was  described, 
implemented  and  applied  to  a  dataset.  Comparison  of  the  reconstructed  3D  shape  to  ver¬ 
ification  data  showed  that  the  reconstruction  was  smoother  and  more  rounded  than  the 
verification  data.  Despite  the  smoothing,  it  is  possible  to  reconstruct  the  grain  with  an 
average  error  of  2.7  micrometers  on  the  location  of  the  grain  boundary. 
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